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Abstract 

We calculate corrections to the fermion propagator and to the Green's functions of all fermion 
bilinear operators of the form \&r\I r , to one-loop in perturbation theory. 

We employ the Wilson/clover action for fermions and the Symanzik improved action for gluons. 

The novel aspect of our calculations is that they are carried out to second order in the lattice 
spacing, 0{a 2 ). Consequently, they have addressed a number of new issues, most notably the 
appearance of loop integrands with strong IR divergences (convergent only beyond 6 dimensions) . 
Such integrands are not present in 0(a 1 ) improvement calculations; there, IR divergent terms are 
seen to have the same structure as in the O(a ) case, by virtue of parity under integration, and 
they can thus be handled by well-known techniques. We explain how to correctly extract the full 
0{a 2 ) dependence; in fact, our method is generalizable to any order in a. 

The 0(a 2 ) corrections to the quark propagator and Green's functions computed in this paper 
are useful to improve the nonperturbative RI-MOM determination of renormalization constants 
for quark bilinear operators. 

Our results depend on a large number of parameters: coupling constant, number of colors, lattice 
spacing, external momentum, clover parameter, Symanzik coefficients, gauge parameter. To make 
these results most easily accessible to the reader, we have included them in the distribution package 
of this paper, as an ASCII file named: Oa2results.m ; the file is best perused as Mathematica input. 

PACS numbers: 11.15.Ha, 12.38.Gc, ll.lO.Gh, 12.38.Bx 

Keywords: Lattice QCD, Lattice perturbation theory, Fermion propagator, Fermion bilinears, Improved 
actions 
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I. INTRODUCTION 



A major issue facing Lattice Gauge Theory, since its early days, has been the reduction 
of effects which are due to the finite lattice spacing a, in order to better approach the elusive 
continuum limit. A systematic framework to address this issue is Symanzik's program 
in which the regularized action is improved through a judicious inclusion of irrelevant 
operators with increasing dimensionality. Thus far, most efforts have been directed towards 
0(a}) improvement; this is automatic in some cases (i.e. requires no tuning of parameters), 
by symmetry considerations alone. Such is the case, for example, of the twisted mass 
formulation of QCD at maximal twist, where certain observables are 0(a l ) improved, 
as a consequence of symmetries of the fermion action: Setting the maximal twist requires 
the tuning of only a single parameter in the action, i.e. the critical quark mass, and no 
further improvement of the operators is required. 

In other cases, such as with the clover fermion action, 0(a) corrections must be also imple- 
mented on individual operators; such corrections take the form of an additional, finite (non 
UV- divergent) renormalization or an admixture of appropriate higher dimensional opera- 
tors. Determining the values of the renormalization functions or mixing coefficients requires 
an evaluation of appropriate Green's functions, as dictated by the choice of renormalization 
scheme; these Green's functions can be evaluated perturbatively or nonperturbatively. 

As regards the perturbative evaluation of Green's functions for the "ultralocal" fermion 
bilinear operators = ^/\ a T^/ (T denotes all possible distinct products of Dirac matrices, 
and A a is a flavor symmetry generator) and the related fermion propagator, the following 
types of calculations have appeared thus far in the literature: (i) One-loop calculations 
to O(a , In a) have been performed in the past several years for a wide variety^of^ actions, 
ranging from Wilson fermions / gluons to overlap fermions and Symanzik gluons 



9. 



L 1 1 01] (ii) There exist one-loop computations of 0(a 1 ) corrections, with an arbitrary fermion 
mass 



(Hi) The first two-loop calculations of Green's functions for 0„ were completed 
recently, to O(a ), for Wilson/clover /twisted-mass fermions and Wilson gluons 12, Ll3|] . (iv) 



A number of Q(a°) results have also been obtained by means of stochastic perturbation 







Q,ll6|. 



theory 

One-loop computations of 0(a 2 ) corrections did not exist to date; indeed they present 
some novel difficulties, as compared to 0(a}). Extending O(a ) calculations up to O^a 1 ) 



does not bring in any novel types of singularities: For instance, terms which were convergent 
to O(a ) may now develop at worst an infrared (IR) logarithmic singularity in 4 dimensions 
and the way to treat such singularities is well known; also, in most of the cases, e.g. for 
m = 0, terms which were already IR divergent to O(a ) will not contribute to 0(a l ), by 
parity of loop integration. On the contrary, many IR singularities encountered at 0(a 2 ) 
would persist even up to 6 dimensions, making their extraction more delicate. In addition 
to that, there appear Lorentz non-invariant contributions in 0(a 2 ) terms, such as Ylfj,Pfi/P 2 
(p: external momentum). 

In this paper we present a one-loop perturbative calculation, to 0(a 2 ), of the quantum 
corrections to the fermion propagator and to the complete basis of local fermion bilinear 
currents \i/r\I/, using massless fermions described by the Wilson/clover action. We use a 
3-parameter family of Symanzik improved gluon actions, comprising all cases which are in 



common use (Wilson, tree-level Symanzik, Iwasaki 17|, DBW2 [3], Liischer-Weisz 19l. |20|). 
All calculations have been performed for generic values of the gauge parameter. Also, by 
virtue of working in a massless scheme, all of our results are applicable to other ultralocal 



fermion actions as well, such as the twisted mass or Osterwalder-Seiler action 2jJ. Our 
results can be used to construct 0(a 2 ) improved definitions of the fermion bilinears. In 



particular, they will be used in Ref. [22j to improve the nonperturbative determinations, 



with the RI-MOM method [23[, of renormalization constants of bilinear quark operators. 

This paper is organized as follows: Section II is an outline of our calculational procedure; 
Section III describes in detail the evaluation of a prototype IR divergent integral; Sections 
IV and V present the corrections to the propagator and to fermion bilinears, respectively; 
Section VI contains a discussion and concluding remarks. Appendix A contains a basis of 
the divergent integrals which appear in the calculation, evaluated to the required order in 
a. 
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II. DESCRIPTION OF THE CALCULATION 



Our calculation makes use of the clover (SW) action for fermions; for Nf flavor species 
this action reads, in standard notation, 

S F = ^2^2(4r + m f )il}f(x)rl}f(x) 
f x 



ip f (x) (r - 7jU ) U x , x+fJ ,ip f (x + fi)+ip f (x + n) (r + 7/ J U x+ ^ x ip f {x) 

The Wilson parameter r is set to r = 1; / is a flavor index; = [7^, 7„]/2; the clover 
coefficient csw is kept as a free parameter throughout. Powers of the lattice spacing a have 
been omitted and may be directly reinserted by dimensional counting. The tensor is a 
lattice representation of the gluon field tensor, defined through 

Ffiv = g (QflU ~ Ql>[i) (2) 

where Q^ u is the sum of the plaquette loops 

Qflf U x ,x+IJ,Ux+fl, X+/J,+uU x +fl+U, X + uU X -\-U,X ~\- Ux y x + uUx + U,X + U — /J,Ux + U — fl,X — llUx — fj -y x 

U Xl x—^Ux—fjL, x—^—yUx—^—v, x—vUx—u, x ~\~ U Xy x— vU x — V) x— V+jJJx— U+fJ,, x+fiUx+^i, x (3) 

We perform our calculation for mass independent renormalization schemes, so that mj = 
0; this simplifies the algebraic expressions, but at the same time requires special treatment 
when it comes to IR singularities. By taking m/ = 0, our calculation and results are identical 
also for the twisted mass action and the Osterwalder-Seiler action in the chiral limit (in the 
so called twisted mass basis). 

For gluons we employ the Symanzik improved action, involving Wilson loops with 4 and 
6 links 1 , which is given by the relation 

2 r ^ ^ 

c ^ ReTr{l - f/piaqucttc} + ci Re Tr {1 - £/ rcct angic} 

plaquette rectangle 

+ c 2 Re Tr { 1 - £/ chair } + c 3 ^ Re Tr { 1 - £/ para iieiogram } (4) 

chair parallelogram 



$g — -9 



1 1 x 1 plaquette, 1x2 rectangle, 1x2 chair (bent rectangle), and lxlxl parallelogram wrapped around 
an elementary 3-d cube. 
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The coefficients q can in principle be chosen arbitrarily, subject to the following normaliza- 
tion condition, which ensures the correct classical continuum limit of the action 

c + 8ci + 16c 2 + 8c 3 = 1 (5) 

Some popular choices of values for q used in numerical simulations will be considered in 
this work, and are itemized in Tabled! they are normally tuned in a way as to ensure 0(a 2 ) 
improvement in the pure gluon sector. Our one-loop Feynman diagrams do not involve 
pure gluon vertices, and the gluon propagator depends only on three combinations of the 
Symanzik parameters: Co = Co + 8ci + I6C2 + 8C3 (= 1), G\ = c-i + C3, C2 = c\ — C2 — C3 ; 
therefore, with no loss of generality all these sets of values have C2 = 0. 

For the algebraic operations involved in evaluating the Feynman diagrams relevant to 
this calculation, we make use of our symbolic package in Mathematica. Next, we briefly 
describe the required steps: 

• Algebraic manipulations: 

The first step in evaluating each diagram is the contraction among vertices, which is per- 
formed automatically once the algebraic expression of the vertices and the topology ( "inci- 
dence matrix" ) of the diagram are specified. The outcome of the contraction is a preliminary 
expression for the diagram; there follow simplifications of the color dependence, Dirac ma- 
trices and tensor structures. We also fully exploit symmetries of the theory (periodicity, 
reflection, conjugation, hypercubic, etc.) to limit the proliferation of the algebraic expres- 
sions. 

• Dependence on external momentum: 

Even though one-loop computations are normally a straighforward procedure, extending 
to 0(a 2 ) introduces several complications, especially when isolating logarithms and Lorentz 
non-invariant terms. As a first task we want to reduce the number of infrared divergent inte- 
grals to a minimal set. To do this, we use two kinds of subtractions among the propagators, 



6 



using the simple equalities 



1 



4E,sm 4 (g,/2)-4(E M sm 2 (g,/2) 



, , + -T^K (6) 

q l q A q z q z 

D(q) = D plaq (q) + j/J(g) - D plm (q)} 

= D plaq (q) + D plaq (q){D; l 1 aq (q) - D-\q))D{q) (7) 

where q stands for k or k + ap, and & (p) is the loop (external) momentum. The denominator 
of the fermion propagator, q 2 , is defined as 

f = ^sin%)+(m / + ^ 2 ) 2 , g 2 = 4^sin 2 (|) (8) 



Q 



For the present work, one sets rrif = and r = 1, as used in Eq. (Q; D is the 4x4 
Symanzik gluon propagator; the expression for the matrix (DZ^q) — D^ 1 (q)), which is 
0(q 4 ), is independent of the gauge parameter, A, and it can be easily obtained in closed 
form. Moreover, we have 

(*W?))„„ = S f - (1 - A)^ 0) 

Terms in curly brackets of Eqs. ([6]) and ((7|) are less IR divergent than their unsubtracted 
counterparts, by two powers in the momentum. These subtractions are performed iteratively 
until all primitively divergent integrals (initially depending on the fermion and the Symanzik 
propagator) are expressed in terms of the Wilson gluon propagator. 

Having reduced the number of distinct divergent integrals down to a minimum, the most 
laborious task is the computation of these integrals, which is performed in a noninteger num- 
ber of dimensions D > 4. Ultraviolet divergences are explicitly isolated a la Zimmermann 
and evaluated as in the continuum. The remainders are D-dimensional, parameter-free, zero 
external momentum lattice integrals which can be recast in terms of Bessel functions, and 
finally expressed as sums of a pole part plus numerical constants. We analytically evaluate 
an extensive basis of superficially divergent loop integrals, listed in Eqs. ( 1A1I) - ( 1A10I) of 



Appendix El a few of these were calculated in Ref. |2J]. The integrals of Eqs. (TAT]) . (lA"2l . 
(1A3I) . are the most demanding ones in the list; they must be evaluated to two further or- 
ders in a, beyond the order at which an IR divergence initially sets in. As a consequence, 



their evaluation requires going to D > 6 dimensions. Fortunately, they are a sufficient basis 
for all massless integrals which can appear in any 0(a 2 ) one-loop calculation; that is, any 
such computation can be recast in terms of (lAip . (1A2I) . (IA30 . plus other integrals which are 
more readily handled. A correct way to evaluate flAip . (IA2I) . (IA3[) has not been presented 
previously in the literature, despite their central role in 0(a 2 ) calculations, and this has 
prevented one-loop computations to 0(a 2 ) thus far. The calculation of such an integral is 
given in detail in the next section. 

Terms which are IR convergent can be treated by Taylor expansion in ap to the desired 
order. Alternatively, the extraction of the ap dependence may be performed using iteratively 
subtractions of the form 

f(k + ap) = f(k) + [f(k + ap) - /(*)] (10) 
This leads to exact relations such as the following ones 

1 1 E A1 sin ( 2fc M + a ?V) Sill ( a ?V) 



k + ap ^ k + ap k 2 



E M sin(A; M + sin(^) (k 2 + k + 7p 
k+ap k 2 



(11) 



1 = J_ _ 4E / ,sin(fc M + ^)sin(^) 
k + ap k 2 k + ap k 2 

In these relations the exact ap dependence of the remainders is under full control; this type 
of subtraction is especially useful when applied to the Symanzik propagator. 

• Numerical integration: 

The required numerical integrations of the algebraic expressions for the loop integrands (a 
total of ~ 40,000 terms) are performed by highly optimized Fortran programs; these are 
generated by our Mathematica 'integrator' routine. Each integral is expressed as a sum 
over the discrete Brillouin zone of finite lattices, with varying size L (4 4 < L A < 128 4 ), and 
evaluated for all values of the Symanzik coefficients listed in Table [I] (corresponding to the 
Plaquette, Symanzik, Iwasaki, TILW and DBW2 action). 
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• Extrapolation: 

The last part of the evaluation is the extrapolation of the numerical results to infinite 
lattice size. This procedure entails a systematic error, which is reliably estimated, using a 
sophisticated inference technique; for one-loop quantities we expect a fractional error smaller 
than 1(T 7 . 



III. EVALUATION OF A PRIMITIVELY DIVERGENT INTEGRAL 

Divergent integrals which appear in calculations up to 0(a l ) may be evaluated using the 



standard procedure of Kawai et al. [25j, in which one subtracts and adds to the original 
integrand its naive Taylor expansion, to the appropriate order with respect to a, in D — > 4 + 
dimensions: The subtracted integrand, being UV convergent, is calculated in the continuum 



limit a — > 0, using the methods of Ref. 26(], while the Taylor expansion terms are recast in 
terms of Bessel functions and are evaluated in the limit e — > (e = (4 — D)/2). 

In contrast to the above, some of the integrals in the present work, given that they must 
be evaluated to 0(a 2 ), have Taylor expansions which remain IR divergent all the way up 
to D < 6 dimensions. A related difficulty regards Kawai's procedure: Subtracting from the 
original integral its Taylor expansion in .D-dimensions to the appropriate order, the UV- 
convergent subtracted expression at which one arrives can no longer be evaluated in the 
continuum limit by naively setting a — > 0, because there will be 0(a 2 ) corrections which 
must not be neglected. These novel difficulties plague integrals IA11 IA21 \A3\ of Appendix [A] 
Using a combination of momentum shifts, integration by parts and trigonometric identities, 
one may express IA2I and IA3I in terms of IA1I and other less divergent integrals. Thus, it 
suffices to address the evaluation of IA1I 



MP) = I 7^,— l^nr (13) 
W k 2 k + ap 



This is a prototype case of an integral which is IR divergent in D < 6 dimensions; in fact, all 
other integrals encountered in the present calculation may be expressed in terms of Al(p) 
plus other integrals which are IR convergent at D > 4 (and are thus amenable to a more 
standard treatment). 
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First we split the original integrand / into two parts 

1= : 2~^2 =h + h (14) 
k 2 k + ap 

where I2 is obtained from J by a series expansion, with respect to the arguments of all 
trigonometric functions, to subleading order; Ii is simply the remainder I — 12 

k 2 -^-k 2 k 4 (^k 2 -k 2 ^J k A ((k + ap) 2 - k~+a~p\ 

^1 = 2 "I 2 "I 2 

k 2 k 2 k + ap 12 (A; 2 ) 2 k 2 k + ap 12 (A; 2 ) 2 (k + ap) 2 k + ap 

(k + ap) 2 - - fc + ^> 2 ( fc + a P) 4 (i k + ap) 2 ~ k + 7p 2 ) 

+ " " 12 ^-^ 2 + " 2 ^-_ 2 7 (15) 

k 2 (k + ap) 2 k + ap 12 A; 2 ((k + ap) 2 ) k + ap 



1 

+ 



(A; + ap) 4 



+ 



12 A; 2 ((Jfc + ap) 2 ) 2 12 (A; 2 ) 2 (Jfc + ap) 2 



(16) 



A; 2 (A; + ap) 2 

(g 4 = Yl^Qfi)- 1% * s f ree °f trigonometric functions, while I\ is naively Taylor expandable to 
0{a 2 ); its integral equals 

^ d A k 



-h = 0.004210419649(1) + a 2 p 2 0.0002770631001(3) + 0(a 4 , a 4 In a) (17) 

-7T ( 2 ^) 4 

The errors appearing in the above equation come from extrapolations to infinite lattice size. 

To evaluate the integral of I2 we split the hypercubic integration region into a sphere of 
arbitrary radius /i about the origin (fi < 7r) plus the rest 



+ / "/ (18) 

The integral outside the sphere is free of IR divergences and is thus Taylor expandable to 
any order, giving 2 (for /1 = 3.14155) 

(/" - / ) (0I J 2 = 6.42919(3) 10- 3 + a 2 p 2 6.2034(1) lO" 5 + 0(a A ) (19) 

We are now left with the integral of I2 over a sphere. The most infrared divergent part of 
J 2 is 1/(A; 2 {k + ap) 2 ), with IR degree of divergence -4, and can be integrated exactly, giving 

"** 1 1 'l-ln(^)) (20) 



2 Due to its peculiar domain, this integral has been evaluated by a Monte Carlo routine, rather than as a 
sum over lattice points. The errors in Eq. (|19p are thus Monte Carlo errors. 
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The remaining two terms comprising I2 have IR degree of divergence -2, thus their calculation 
to 0(a 2 ) can be performed in D-dimensions, with D slightly greater that 4. Let us illustrate 
the procedure with one of these terms: k A /((k 2 ) 2 (k + ap) 4 ). By appropriate substitutions 
of 

1 _ 1 -2(k-p) -f 
~ F 4 



(k + p) 2 
we split this term as follows 

A- 1 



k 2 (k+p) 2 



(p = ap) 



(21) 



{k 2 ) (k + p) 2 [(k 2 ) 



k A k 4 (-2(k ■ p) - p 2 ) Ak\k-p) 2 

t 77m y 



(k 2 



+ 



k A (4 (A; • p)p 2 + (p 
(A; 2 ) 4 (k + p) 2 



+ 



{k 2 f 

Ak\k ■ pf {-2{k ■ p) - p 2 ) 
(A; 2 ) 5 (k + p) 2 



(22) 



The part in square brackets is polynomial in a and can be integrated easily, using D- 
dimensional spherical coordinates. The remaining part is UV-convergent; thus the integra- 
tion domain can now be recast in the form 



\k\</j, 



(23) 



fe|<oo J ii<\k\<oo 



The integral over the whole space can be performed using the methods of Ref. 26(, whereas 
the integral outside the sphere of radius fi is 0(a 3 ) and may be safely dropped. The same 
procedure is applied to the last term of 12- Adding the contributions from all the steps 
described above, we check that the result is independent of \i. 



IV. CORRECTION TO THE FERMION PROPAGATOR 

The fermion propagator is the most common example of an off-shell quantity suffering 
from O(a) effects. Capitani et al. j(| have calculated the first order terms in the lattice spac- 
ing for massive fermions. We carried out this calculation beyond the first order correction, 
taking into account all terms up to 0(a 2 ). Our results, to O^a 1 ), are in perfect agreement 
with those of Ref. 6|. The clover coefficient c^w has been considered to be a free parameter 
and our results are given as a polynomial of csw- Moreover, the dependence on the number 
of colors N, the coupling constant g and the gauge fixing parameter A, is shown explicitly. 
The Symanzik coefficients, q, appear in a nontrivial way in the propagator and, thus, we 
tabulate these results for different choices of c,-. 
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The one-loop Feynman diagrams that enter our 2-point Green's function calculation, are 
illustrated in Fig. 1. 




1 2 
Fig. 1: One-loop diagrams contributing to the fermion propagator. Wavy (solid) 
lines represent gluons (fermions). 



Next, we provide the total expression for the inverse fermion propagator S^ 1 as a function 
of g, N, csw; A. Here we should point out that for dimensional reasons, there is a global 
prefactor 1/a multiplying our expressions for the inverse propagator, and thus, the 0(a 2 ) 
correction is achieved by considering all terms up to 0(a 3 p 3 ). 



' T 6 



if £ (0,1) - 4.79200956(5) A + e (0 ' 2) c sw + e m 4 W + A ln(aV 



2 ~2 

- ap g 



- 3.86388443(2) A + e (1 ' 2) c sw + £ (1,3) 4 W - - (3 - 2 A - 3 c sw ) ln(aV) 



- %a 2 f g 2 £ (2 ' 1} + 0.507001567(9) A + s {2 > 2) c sw + e^'°' Cg W 



,(2,3) J2 



+ 



101 _ 11 _A 
120 30 2 6 



ln(a 2 p 2> 



- ia 2 p 2 jg 2 e (2A) + 1-51604667(9) A + e< 2,5) c sw + ^ {2 ' 6) c 2 w 



, 59 ci C 2 1 /3, 2 

h 1 1 -A + c sw + c 2 

240 2 60 4 I 2 bW 



sw 



2„2\ 



In (a p 



- jo p z g 



i _ 5k _ 

30 10 48 J 



(24) 



where t? 2 = (? 2 C F /(167r 2 ), C F = (iV 2 - l)/(2iV), C 2 = c x - c 2 - c 3 , f = E M 7^ and 
the specific values A = 1 (0) correspond to the Feynman (Landau) gauge. The quantities 
e^'fi appearing in our results for S^ 1 are numerical coefficients depending on the Symanzik 
parameters, calculated for each action we have considered and tabulated in Tables HT1 - IVt 
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the index % denotes the power of the lattice spacing a that they multiply. In all Tables, the 
systematic errors in parentheses come from the extrapolation over finite lattice size L — > oo. 

Terms proportional to 1/a have been left out of Eq. (1241) for conciseness; such terms 
represent 0(g 2 ) corrections to the critical value of the fermion mass. 

We observe that the 0(a}) logarithms as well as all terms multiplied by A, are independent 
of the Symanzik coefficients; on the contrary 0(a 2 ) logarithms have a mild dependence on 
the Symanzik parameters. A number of Lorentz non-invariant tensors (Y^^pfii / 3 ) appear in 
O (a 2 ) correction terms, compatibly with hypercubic invariance. Finally, our 0{a}) results 
for the Plaquette action, are in agreement with Eq. (37) of Ref. jf]]. 

To enable cross-checks and comparisons, the per-diagram contributions d\(p), ^(p) are 
presented below. The tadpole diagram 1 of Fig. 1 is free of logarithmic terms and indepen- 
dent of csw; its final expression is 



di(p) 
~9 2 



=(o,i) 



+ 3.050262540200(1) A 



+ ap 



=(i,i) 



+ 1.529131270100(1) A 



+ % a 



2 Ji 



=(2,1) 



- 0.509710423367(1) A 



(25) 



where the numerical values for the Symanzik dependent coefficients if are listed in Ta- 
ble IVB The main contribution to the propagator correction comes from diagram 2, as can be 
seen from the following expression, with provided in Table IVIII The remaining terms 
with coefficients e^'i* are the same as in Eq. (T5] 
d 2 (p) 



~9 2 



ef l) - 7.850272109(6) A + e [V ' 2) c sw + 4 W + A \n(a 2 p 



(0,2) 



,(0,3) 2 



+ ap 



=(i,i) 



5.39301570(2) A + e (1 ' 2) c sw + e 



2J2\ 



+ ia 



2 Ji 



£ ~(2,D 



+ 



., + 1.016711991(9) A + e (2 ' 2 > c sw + e 
11 



(2,3) .2 



(3 — 2 A — 3 csw) m (a p 



sw 



101 
120 



30 



ln(a 2 p 2> 



+ i a 2 p 2 p- e {2 ^ + 1.51604667(9) A + c sw + c 



(2,5) 



r(2,6) -2 



SW 



+ i a 



+( 


f 59 
v 240 + 


2 


c 2 

60 " 




r 3 


c 2 


5 


p 2 


I 80 


10 


48 



1 fi 



4 V 2 



-sw 



2J2\ 



In (a p 



A 



(26) 
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Using our results for the fermion propagator, we can compute the multiplicative renor- 
malization function of the quark field (Z^). 

V. FERMION BILINEAR OPERATORS 

In the context of this work we also study the 0(a 2 ) corrections to Green's functions of 
local fermion operators that have the form ^r^. T corresponds to the following set of 
products of the Dirac matrices 

r = 1, 7 5 , 7/1 , 7 5 7 M , 7 5 cv, = ^[ 7/1 , lv \ (27) 

for the scalar {0 s ), pseudoscalar {O p ), vector {O v ), axial (0 A ) and tensor (0 T ) operator, 
respectively. We restrict ourselves to forward matrix elements (2-point Green's functions, 
zero momentum operator insertions). We also considered the tensor operator T ', corre- 
sponding to T = a^y and checked that the Green's function coincides with that of T ; this 
is a nontrivial check for our calculational procedure. 

The only one-particle irreducible Feynman diagram that enters the calculation of the 
above operators is shown in Fig. 2. 




Fig. 2: One-loop diagram contributing to the bilinear operators. A wavy (solid) 
line represents gluons (fermions). A cross denotes the Dirac matrices 1 (scalar), 7 s 
(pseudoscalar), 7^ (vector), 7 5 7 M (axial), 7 5 <r (UJ , (tensor T) and (tensor T"). 

We show our results for the one-loop corrections to the amputated 2-point Green's func- 
tion of each operator ^/T^/, at momentum p 

A r ( p ) = r#) q)^ p (28) 

Our final results are given as a polynomial of Csw 5 in a general covariant gauge. Since 
their dependence on the Symanzik parameters, q, cannot be written in a closed form, as 
in the case of the quark propagator we will tabulate the numerical coefficients for a variety 
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of choices for q, in order to cover a range of values that are used in both perturbative 
calculations and numerical simulations. 

We begin with the 0(a 2 ) corrected expression for A s (p); including the tree-level term, 
we obtain 



(P) 



l + ~9 2 



4 0,1) + 5.79200956(5) A + ef z> c sw + e£' aj c 2 sw - ln(aV) (3 + A) 

(1=3) 9 



-(0,2) 



-(0,3) 7 



7.2 



+ aip'g 



ef x) - 3.93575928(1) A + 4' 2) c sw + e£ '°> c 2 w + [ ^ + A + ^ c sw 



3 
2 



2 2\ 

a p ) 



, 2 2-2 

+ a p g 



-(2,1) 



2.27358943(5) A + ef 2) c sw + ef*> c 2 w + ( -- A + ^A + £ c sw ) ln(a 2 p : 



-(2,3) 2 



4 4 



*2^2\ 



+ a 2 ^f^ 2 



13 C2 _ A 
24 + T ~ 8 



(29) 



The numerical coefficients £5 , £5^ and e^'^ with their systematic errors are presented in 
Tables IVffll HXl and H respectively. 

One might attempt to use the 0(a) corrections computed above in order to devise an 
improved operator, with suppressed finite-a artifacts; it should be noted, however, that 
improvement by means of local operators, as permitted by Quantum Field Theory, is not 
sufficient to warrant a complete cancellation of 0(a 2 ) terms in Green's functions, since the 



latter contain also terms with non-polynomial momentum dependence, such as Yl u Pt/p 2 



v.? in 



Thus, at best, one can achieve full 0(a 2 ) improvement only on-shell, or approximate im- 
provement near a given reference momentum scale. Such non-polynomial terms are not 
present at 0(a x ). This comment applies also to the remaining operators we examine below. 

Next, we turn to A p (p), where Symanzik dependent coefficients, £p \ are tabulated in 
Table IXII The pseudoscalar operator is free of 0(a l ) terms; moreover, all contributions 
linear in csw vanish 



A'(p) = 7 5 + 7 5 ^ 2 



1 2 2 5-2 

+ a P 19 



-(0,1) 



-(2,1) 



5.79200956(5) \ + e 



(0,2) 2 



. - 0.83810121(5) A + e {2 ' 2) c 2 sw + ( ~ + ~ ) ln{a 2 p 



ln(aV) (3 + A) 



1 A 



2„2> 



+ a 2 ^f^~g 2 



13 C 2 Ai 
24 ~2 ~ ! 



(30) 
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The 0(a 2 ) corrected expressions for A v (p) and A A (p) are more complicated, compared 
to the scalar and pseudoscalar amputated Green's functions, in the sense that momentum 
dependence assumes a variety of functional forms; this fact also introduces several coefficients 
which depend on the Symanzik parameters 



2A 



+ 7„r 



4 0,1) + 4.79200956(5) A + ep 2) c sw + 4°' 3) 4w ~ A ln(aV)] 
e { y' 1} - 0.93575928(1) A + e ( y' 2) c sw + 4' 3) c sw 



+ (-3 + A + 3 c sw ) ln(aV) 



i 2 2-2 



-(2,1) , A 



■-v 4 2 ' 2) csw + 4 2 ' 3) 4w + [~ + WaV) 



+ ^ 7^ 2 f [4 2 ' 4) " 0.8110353(1) A + 4' 5) csw + # 6) 4v 



11 ci C 2 A 5 



+ a 2 j^0 2 [e£ ,7) + 0.2436436(1) A + eg >8) c sw + eg' 9) c| w 



+ 



y p A 

2 ~2 

+ a 



Ci A csw . c 



149 , , ^ 

'120 Cl 30 + 4 + 6 + 



'sw 



ln(aV) 



3 C 2 5 
80 + 10 + 48 . 



+ a —^-g 



1 2^ Ai 

h -C 2 H 

60 5 12 



+ a 



2 



y p 4 



,2\2 3 



101 11 A 

"60" + 15 C2+ 3 
3 C 2 5 



40 



24 



A 



(31) 
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The numerical values of ep for different Symanzik choices are given in Tables Kill - IXVI1 



A A (p) = 7 5 7m + 



2A 



TPPf* -2 

— 9 

ef 1] + 4.79200956(5) A + ef 2) c sw + 4 W - A ln(a 2 p 



-(0,3) 9 



7.2 



e^ 1] - 2.93575928(1) A 



+ £a' 2) c sw + ^a'*' Cs W + * ln(° 2 P 



(1.3) J2 



,2^2\ 



, 2 5 2-2 



, 2 5 2 -2 

+ a 7 7 M P 5 



-(2,1) 



A 



-(2,2) 



-(2,3) 2 



-+e A ' csw + e^ c sw + ^ i2Q , iQ 

4 2 ' 4) - 1.7465235(1) A + e% 5) c sw + 4 2 ' 6) 4w 
109 d C 2 5 X 7 



53 11 



+ ^C 2 ) ln(aV^ 



+ 



~ TT ~ 7^77 + o A + 77: C SW - 



sw 



240 2 60 8 12 4 
- a 2 7 5 /P^ 2 + 1-1146200(1) A + 4 2 ' 8) c sw + 4 2 ' 9) 4w 



ln(a 2 p 2> 



+ 



91 C 2 3 5 Cg W 

ci A c S w — 

120 30 4 6 2 



ln(a 2 p 2 ^ 



2 5 t—iprp ~2 

+ a 7 7 M — 



p- 



. 2 5 ~2 

+ a 1—T9 



80 10 48 . 



1 2^ A 

h -C 2 H 

60 5 12 



, 2 5 M^- 2 
+ a 7 ,,. 2 x 2 g 



(P 



101 11 A 
~6¥ + 15 2+ 3J 

_ — - 9i - — 

~ 40 ~ T ~ 24 



(32) 



Eq. (j3ip and Eq. (1321) have many similar terms, among them the coefficients 



-(0,2) 



(0,2) (0,3) 



— e 



(0,3) 



(33) 



The rest of the coefficients el appear in Tables IXVIII - IXXI 



The remaining Green's functions that we computed are those corresponding to the tensor 
bilinears (T = 7 5 a^, T' = a^), which are the most complicated of all the operators that 
we studied. Clearly, the Green's functions A T (p) and A T '(p), corresponding to T and T", 
coincide numerically, even though this fact is not immediately apparent from their algebraic 
forms. In fact, we computed both A T (p) and A T '(p) in two distinct calculations; their 
numerical coincidence constitutes a rather nontrivial check of our results. For the reader's 
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convenience, we present below both tensor Green's functions. 



A T (p) = 7 5 <V + 7 5 <V 9 Z 4°' 1} + 3.79200956(5) A + e^' c sw + 4' a; 4 W + (1 - A) \n(a z p 



(0,2) 



-(0,3) 2 



2 2^ 



+ GU7 g 



4 1,1} + 3.87151852(5) A + e£' 2) c sw + 4' 3) 4w 
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2 2. 
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2p 2 



.9" 
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17 a 

L 6 3 
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1 C 2 A 
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7 2p 2 P 
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(34) 



The coefficients 4 are tabulated in Tables IXXII - IXXV1 
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A r {p) = <v + g 2 + 3.79200956(5) A + eJ?/ 2) c sw + 4°' 3) c sw + (1 - A) ln(a 2 p 2 ) 

4V } - 3.87151852(5) A + e£; 2) c sw + 4' 3) 4 W 
- (-3 + 2A + c sw ) ln(aV) 



+ a' 
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(35) 



Several coefficients can be written in terms of Et (Eqs. (|36|) - (1551) ). while the rest are 
given in Tables IXXVTl - IXXVIII 
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(36) 
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"T J 
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-(2,5) 



.(2,6) 

Cy/ 



-(2,6) 



(37) 
(38) 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have calculated the fermion propagator S(p) and the Green's functions 
A r (p) for the fermion bilinear operators ^DD, where T stands for any product of Dirac 
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gamma matrices. Our calculations were performed to one loop in lattice perturbation theory, 
using the Wilson/clover fermion action. For gluons we employed a family of Symanzik 
improved actions, parameterized by 3 independent "Symanzik" coefficients; explicit results 
are presented for some of the most commonly used actions in this family: Wilson, Tree-level 
Symanzik, Tadpole improved Liischer-Weisz, Iwasaki and DBW2. 

Our calculations extend, to a rather large family of fermion/gluon actions, results which 
were previously known to O(a ) and 0{a 1 ) (modulo In a). However, the truly novel fea- 
ture in our calculations is that they were performed to second order in the lattice spacing 
a (O (a 2 , a 2 In a)). This fact introduces a number of complications, which are not present 
in lower order results. In a nutshell, the reason for these complications is as follows: The 
extraction of a further power of a from a Feynman diagram strengthens, by one unit, the 
superficial degree of infrared (IR) divergence of the corresponding integrand over loop mo- 
menta. Thus, a priori, in a 0(a l ) calculation, loop integrals would be IR convergent only 
in D > 5 dimensions; however, as can be easily deduced by inspection, the most divergent 
parts of the integrands are odd functions of the loop momenta, and will thus vanish upon 
integration. What is left behind is a less divergent integrand which is IR convergent in 
D > 4, just as in the case of O(a ) calculations, and can thus be treated by standard meth- 



ods, such as those of Ref. 



251 ] . For 0(a 2 ) calculations, on the other hand, integrands are IR 



convergent only at D > 6, and their most divergent parts no longer vanish upon integration; 
a naive application of the procedure of Ref. 25] will fail to produce all 0(a 2 ) contributions. 
The procedure which we propose in this work for handling the above difficulty is in fact 
applicable to any order in a. In brief, it recasts the integrands as a sum of two parts: The 
first part can be exactly evaluated as a function of a, while the second part is naively Taylor 
expandible, as a polynomial to the desired order in a. 

Since the propagator and Green's functions are meant to be used in mass independent 
renormalization schemes, our results have been obtained at vanishing fermionic masses; the 
case of massive fermions (including non-degenerate flavors and twisted mass terms) will 
appear in a forthcoming publication. Nevertheless, even at vanishing masses, our final ex- 
pressions are quite lengthy, since they exhibit a rather nontrivial dependence on the external 
momentum (p), and they are explicit functions of the number of colors (N), gauge parame- 
ter (A), lattice spacing (a), clover coefficient (csw) an d coupling constant (g); furthermore, 
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most numerical coefficients in these expressions depend on the Symanzik parameters of the 
gluon action, and we have tabulated them for the actions we have selected. For convenience, 
we accompany this paper with an electronic document, in the form of a Mathematica input 
file, allowing the reader to recover immediately numerical values for any choice of input 
parameters. 

One possible use of our results is in constructing improved versions of the operators O r , 
with reduced lattice artifacts. In doing so, however, one must bear in mind that, unlike 
the 0{a 1 ) case, corrections to 0(a 2 ) include expressions which are non-polynomial in the 
external momentum and, therefore, cannot be eliminated by introducing admixtures of local 
operators. Full improvement can be achieved at best for on-shell matrix elements only. 

Starting from S(p) and A r (p), it is straightforward to write down the renormalization 
functions Z q (for the quark field) and Zr (for the operators O v ) in any renormalization 
scheme. Z q and Z?, as obtained from S(p) and A r (p), differ from the corresponding expres- 
sions evaluated at O(a ), by lattice artefact, which are functions of (a/i) (//: renormalization 
scale), and vanish as a — > 0. At the nonzero values of a employed in numerical simulations, 
these factors are quite important. Ideally, one would prefer a nonperturbative determination 
of renormalization functions; while this is often possible, several sources of error must be 
dealt with. A very effective way to proceed is through a combination of perturbative and 
nonperturbative results. This procedure is carried out and explained in detail in a follow-up 



work 22]]. Briefly stated, nonperturbative data are "corrected" by the perturbative expres- 
sions for Green's functions, and then extrapolated towards small a. As a first illustration 
of this mixed determination, we show in Fig . 3 nonperturbative data for Z q and Zy, deter- 



mined with the RI-MOM method of Ref . [23j , before and after the perturbative corrections. 
The results are obtained by using the Symanzik tree-level improved gluon action at (3 — 3.9 
and the Nf = 2 twisted mass quark action at maximal twist, with gauge field configurations 
and quark propagators generated by the ETM Collaboration 3 . While up to discretization 
effects Zy is a scale independent quantity, the (continuum) RG dependence of Z q on the 
renormalization scale has been removed from the results shown in Fig. 3 by evolving the 
renormalization constant to a fixed reference scale /xq = 1/a (~ 2 GeV), using for the 



3 We thank the members of the ETM Collaboration for having provided us with the data of Fig. 3 before 
publication. 
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anomalous dimension the 4-loop perturbative expression computed in Ref. 27j. Thus, the 
residual dependence of both Z q (fi = 1/a) and Zy on a 2 p 2 observed in Fig. 3 can be safely 
interpreted, at large momenta, as a pure discretization effect. As illustrated in Fig. 3, the 
corrected data are virtually flat, allowing for a safer small-a extrapolation. 




Fig. 3: Non-perturbative data for Z q (iiQ = 1/a) and Zy, before and after perturba- 
tive corrections. Straight lines are extrapolations to small a. (a 2 p 2 = £^ sin 2 (ap /i )) 

The techniques employed in this work are readily applicable to the study of perturbative 

22 



corrections of other Greens's functions, to any desired order in a. Examples are matrix 
elements of 4-fermion operators appearing in effective weak Hamiltonians, and higher di- 
mension twist-2 fermion bilinears involved in generalized parton distributions. We will be 
addressing these issues in forthcoming publications. 
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APPENDIX A: A BASIS OF DIVERGENT INTEGRALS 



The most difficult part of this calculation that requires careful attention is the extraction 
of the dependence on the external momentum p and the lattice spacing a from the divergent 
terms. The singularities are isolated using the procedure explained in section [III an d here we 
present the list of primitively divergent integrals that appeared in our algebraic expressions. 

In the following integrals we define 

k, = 2 sin(^) 

P = 4£sin 2 (^) 

O 

kf, = sin(fc M ) 

In addition, ( ) 5 means sum over inequivalent permutations. No summation over the indices 
/1, v, p, a is implied, unless otherwise stated. 
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Action 




c o 




C3 


Plamipfvfcp 




1.0 


o 


o 


Symanzik 




1 6666667 


-0 083333 

KJ * KJ^JKJKJ KJ 


o 


TILW Ben = . 


S.60 


2 3168064 


-0 151791 


-0 0128098 


TILW /3cn = , 


3.45 


2 3460240 


-0 154846 


-0 0134070 

V7 • KJ .X. t 7 i V t 1 \J 


TILW /3cn = , 


3.30 


2 3869776 


-0 159128 

KJ • ±. KJ kJ A. i-IKJ 


-0 0142442 


TILW Bm — , 


3.20 


2 4127840 


-0 161827 

KJ • J. V 7 1 \J — 1 


-0 0147710 

V7 • W It: 1 I J.KJ 


TILW, [3c = , 


3.10 


2.4465400 


-0.165353 


-0.0154645 


TILW, /5c = . 


3.00 


2.4891712 


-0.169805 


-0.0163414 


Iwasaki 




3.648 


-0.331 





DBW2 




12.2688 


-1.4086 






TABLE I: Input parameters cq, ci, C3. 



Action e (0 ' 1} £ (°' 2 ) £ (°' 3 ) 

Plaquette 16.6444139(2) -2.24886853(7) -1.39726711(7) 
Symanzik 13.02327272(7) -2.01542504(4) -1.24220271(2) 
TILW (8.45) 10.82273528(9) -1.84838009(3) -1.13513794(1) 
TILW (8.00) 10.45668970(6) -1.81821854(5) -1.11582732(3) 
Iwasaki 8.1165665(2) -1.60101088(7) -0.97320689(3) 
DBW2 2.9154231(2) -0.96082198(5) -0.56869876(4) 



TABLE II: The coefficients e^ ^ (Eq. (f24]) ) for different actions. 



29 



Action 

Plaquette 12.8269254(2) -5.20234231(6) -0.08172763(4) 

Symanzik 10.69642966(8) -4.7529781(1) -0.075931174(1) 

TILW (8.45) 9.2865455(1) -4.4186677(2) -0.07160078(1) 

TILW (8.00) 9.0430829(2) -4.35681290(3) -0.070688697(3) 

Iwasaki 7.40724287(1) -3.88883584(9) -0.061025650(8) 

DBW2 3.0835163(2) -2.2646221(1) -0.03366740(1) 



TABLE III: The coefficients (Eq. dMD) for different actions. 



Action 


£ (2,D 


£ (2,2) 


e (2,3) 


Plaquette 


-4.74536466(2) 0.02028705(5) 


0.10348577(3) 


Symanzik 


-4.2478783(2) 


0.05136635(6) 


0.07865292(7) 


TILW (8.45) 


-3.8139475(2) 


0.05751390(9) 


0.06651692(3) 


TILW (8.00) 


-3.7342556(1) 


0.05830392(9) 


0.06444077(4) 


Iwasaki 


-3.2018047(1) 


0.08249970(7) 


0.04192446(4) 


DBW2 


-0.8678072(2) 


0.1024452(2) 


-0.00343999(2) 



TABLE IV: The coefficients e^ 2 ' 1 ) - (Eq. ([21]) ) for different actions. 

Action e {2A) e( 2 < 6 ) 

Plaquette -1.5048070(1) 0.70358496(5) 0.534320852(7) 

Symanzik -1.14716212(5) 0.65343092(3) 0.49783419(2) 

TILW (8.45) -0.92583451(6) 0.62061757(5) 0.467966296(9) 

TILW (8.00) -0.8875297(1) 0.61441084(7) 0.462237852(9) 

Iwasaki -0.6202244(1) 0.55587473(6) 0.41846440(4) 

DBW2 -0.3202477(5) 0.34886590(2) 0.23968038(4) 



TABLE V: The coefficients e (2 ' 4) - e (2 ' 6) (Eq. ([21]) ) for different actions. 
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Action ef> l) 

Plaquette 9.174787621(1) 4.5873938103(5) -1.5291312701(2) 

Symanzik 7.071174701(5) 3.535587351(2) -1.1785291169(8) 

TILW (8.45) 5.86097856(2) 2.930489282(8) -0.976829761(3) 

TILW (8.00) 5.663791993(4) 2.831895997(2) -0.9439653322(7) 

Iwasaki 4.423664730(5) 2.211832365(2) -0.7372774550(8) 

DBW2 1.86908767(4) 0.93454384(2) -0.311514612(6) 

TABLE VI: The coefficients e[ 0,l) (Eq. ([25])) for different actions. 

Action ^ §P ^ 

Plaquette 7.4696262(2) 8.2395316(2) -3.21623339(2) 

Symanzik 5.95209802(7) 7.16084231(8) -3.0693492(2) 

TILW (8.45) 4.96175672(9) 6.3560562(1) -2.8371177(2) 

TILW (8.00) 4.79289770(6) 6.2111869(2) -2.7902902(1) 

Iwasaki 3.6929018(2) 5.19541051(1) -2.4645273(1) 

DBW2 1.0463355(2) 2.1489724(2) -0.5562925(2) 



TABLE VII: The coefficients e£ 0,<) (Eq. ([25])) for different actions. 



Action 


_(o,i) _(0,2) 
b S L s 


_(0,3) 
£ S 


Plaquette 


0.30799634(6) 9.9867847(2) 


0.01688643(6) 


Symanzik 


0.58345905(5) 8.8507071(1) 


-0.12521126(5) 


TILW (8.45) 


0.7049818(1) 8.0538938(2) 


-0.20881716(3) 


TILW (8.00) 


0.7195566(1) 7.9115477(2) 


-0.22196498(3) 


Iwasaki 


0.74092360(2) 6.9016820(2) 


-0.29335071(4) 


DBW2 


-0.0094234(5) 4.0385802(2) 


-0.35869680(4) 



TABLE VIII: The coefficients 4°^ ( E 1- ®) for different actions 
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Action 


e (i,i) 


-(1,2) 
E S 


-(1,3) 


Plaquette 


0.6586287(1) 


-4.20298580(6) 


-1.286053869(4) 


Symanzik 


0.33939970(4) 


-3.76353718(6) 


-1.150059945(4) 


TILW (8.45) 


0.1463203(2) 


-3.42960982(2) 


-1.054472092(1) 


TILW (8.00) 


0.1155729(2) 


-3.36704753(5) 


-1.037165442(1) 


Iwasaki 


-0.05097214(7) 


-2.88571027(1) 


-0.909503374(3) 


DBW2 


-0.1248521(3) 


-1.15247167(2) 


-0.53943631(1) 



TABLE IX: The coefficients 


) (Eq. m: 


) for different actions 


Action ^s'^ 


_(2,2) 
E S 


-(2,3) 
£ S 



Plaquette 2.60041308(7) -4.15080331(7) 0.17641091(2) 

Symanzik 2.3547298(2) -3.85277871(9) 0.196461884(5) 

TILW (8.45) 2.1881285(8) -3.6249313(5) 0.21113016(1) 

TILW (8.00) 2.1605653(8) -3.58171175(4) 0.21385016(2) 

Iwasaki 2.02123300(8) -3.23459547(4) 0.234502732(7) 

DBW2 2.3731619(3) -1.9332087(1) 0.2953480(3) 

TABLE X: The coefficients e ( g' l) (Eq. ([25])) for different actions. 



Action sf^ ef 2) 

Plaquette 9.95102761(8) 3.43328275(3) 0.84419938(7) -0.25823485(3) 
Symanzik 8.7100837(1) 2.98705498(3) 0.70640549(6) -0.27556247(3) 
TILW (8.45) 7.84510495(6) 2.67986902(3) 0.65030355(6) -0.28812231(2) 
TILW (8.00) 7.6896423(1) 2.62578350(2) 0.64432843(6) -0.29027771(3) 
Iwasaki 6.55611308(7) 2.25383382(3) 0.66990790(5) -0.30221183(3) 
DBW2 2.9781769(6) 1.24882665(4) 1.5569125(1) -0.3362271(2) 

TABLE XI: The coefficients and e { p l) (Eq. (f50|>) for different actions. 
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Action Ey 

Plaquette 3.97338480(2) 

Symanzik 3.57961385(3) 

TILW (8.45) 3.32483844(4) 

TILW (8.00) 3.28098129(5) 

Iwasaki 2.98283189(2) 

DBW2 2.25812410(4) 



-(0,2) (0,3) 
t-V V 

-2.49669620(4) 0.85409908(1) 

-2.21267683(2) 0.77806655(1) 

-2.01347343(2) 0.72217154(1) 

-1.97788691(3) 0.71193712(2) 

-1.72542048(4) 0.63679613(2) 

-1.00964505(3) 0.40188086(2) 



TABLE XII: The coefficients e { y l) (Eq. ([51])) for different actions. 



Action 


E v 


-(1,2) 

£ v 


-(1,3) 

E v 


Plaquette 


2.7109817(1) 


-1.84813992(2) 


-0.39052850(2) 


Symanzik 


2.09743725(3) 


-1.51877201(8) 


-0.385127257(2) 


TILW (8.45) 


1.64290440(2) 


-1.2579161(2) 


-0.37793187(2) 


TILW (8.00) 


1.55841933(4) 


-1.20780827(3) 


-0.3761606511(4) 


Iwasaki 


0.9074321(1) 


-0.80352187(4) 


-0.356005234(3) 


DBW2 


-1.4498098(4) 


0.8826550(3) 


-0.264655885(7) 



TABLE XIII: The coefficients e { y l) (Eq. ([51])) for different actions. 



Action 


_(2.1) 

£ v 


-(2,2) 
E V 


_(2,3) 
E V 


Plaquette 


1.5541024(2) 


0.32907377(4) 


-0.0060202576(6) 


Symanzik 


1.6762868(2) 


0.22601986(5) 


0.02822949(2) 


TILW (8.45) 


1.63378530(3) 


0.16772628(3) 


0.04300929(5) 


TILW (8.00) 


1.6190247(1) 


0.15805313(3) 


0.04550457(5) 


Iwasaki 


1.4573118(1) 


0.0858961(2) 


0.07934994(2) 


DBW2 


-1.1604825(4) 


-0.0504803(3) 


0.13992474(3) 



TABLE XIV: The coefficients e { y' 1] - e ( y ,3) (Eq. (j5JJ|) for different actions. 
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Action 


.(2,4) (2,5) 
V 


-(2,6) 
E V 


Plaquette 


0.2500659(2) 0.8859920(1) 


-0.300364436(2) 


Symanzik 


0.0214112(1) 0.8342659(2) 


-0.28736163(1) 


TILW (8.45) 


-0.1100958(1) 0.791026749(4) 


-0.27444757(3) 


TILW (8.00) 


-0.1318272(2) 0.78255494(3) 


-0.27181092(1) 


Iwasaki 


-0.2668492(1) 0.712786719(6) 


-0.25078366(2) 


DBW2 


-0.1528741(6) 0.42190739(4) 


-0.13978037(7) 



TABLE XV: The coefficients e { y A) - e ( y ,6) (Eq. §2)) for different actions. 



ActiOn Ey'^ Ey'^ £y'^ 

Plaquette 1.27887765(9) 0.27776135(2) -0.35475044(2) 

Symanzik 1.03773908(9) 0.28969451(4) -0.302816648(5) 

TILW (8.45) 0.89400856(7) 0.2930984(2) -0.25886703(3) 

TILW (8.00) 0.87034685(8) 0.29343883(9) -0.25031691(5) 

Iwasaki 0.76263373(2) 0.29755270(5) -0.184270928(8) 

DBW2 1.7371355(5) 0.2960594(1) 0.10831780(4) 



TABLE XVI: The coefficients e { y J) - e ( y ,9) (Eq. (JED) for different actions. 



Action 


-(0,1) 
e A 




-(1,2) 
e A 


Jl,3) 
e A 


Plaquette 


-0.84813073(8) 


1.34274645(8) 


-1.71809242(4) 


0.130176166(7) 


Symanzik 


-0.48369852(8) 


0.92541220(1) 


-1.54604828(4) 


0.128375752(4) 


TILW (8.45) 


-0.2452231(1) 


0.64518173(2) 


-1.42093097(4) 


0.125977289(1) 


TILW (8.00) 


-0.20406156(8) 


0.59652190(3) 


-1.39831769(4) 


0.125386884(3) 


Iwasaki 


0.0752372(1) 


0.2684958(1) 


-1.238019617(7) 0.1186684108(9) 


DBW2 


0.7643240(1) 


-0.56650487(5) 


-0.75581589(7) 


0.088218628(3) 



TABLE XVII: The coefficients and e^ i] (Eq. fl32])) for different actions. 
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Action 


-(2,1) 
e A 


-(2,2) 
e A 


-(2,3) 
E A 


Plaquette 


0.3879068(1) 


1.85116980(8) 


-0.093094486(8) 


Symanzik 


0.29616583(7) 


1.7629637(2) 


-0.11136345(3) 


TILW (8.45) 


0.2483248(2) 


1.65782471(1) 


-0.118658311(8) 


TILW (8.00) 


0.2378781(2) 


1.63840648(4) 


-0.11982438(1) 


Iwasaki 


0.05917686(4) 


1.5707047(2) 


-0.13932655(1) 


DBW2 


-2.2341918(4) 


1.22932319(6) 


-0.17119304(8) 



TABLE XVIII: The coefficients - ej' 3) (Eq. $Z§) for different actions. 



Action ^a' 4 '' £ a'^ £ a' 6 ^ 

Plaquette 1.6350438(1) -1.59945524(6) 0.333900263(8) 

Symanzik 1.3008790(1) -1.48761993(3) 0.314172576(5) 

TILW (8.45) 1.0461303(2) -1.39361896(4) 0.297787700(5) 

TILW (8.00) 0.9998744(1) -1.37577372(4) 0.29455820(2) 

Iwasaki 0.6845753(1) -1.24800562(3) 0.26827353(2) 

DBW2 0.0967251(2) -0.735419342(8) 0.14738921(4) 



TABLE XIX: The coefficients ej* 4) - ef 6) (Eq. ([32])) for different actions. 



Action e^' 7 -* e^' 8 ' 1 e^' 9 ^ 

Plaquette 0.41758917(4) 0.395847810(9) 0.31972188(2) 

Symanzik 0.596637529(2) 0.33473715(4) 0.27870681(2) 

TILW (8.45) 0.73021636(8) 0.29171961(4) 0.24115534(2) 

TILW (8.00) 0.75716237(6) 0.28297665(3) 0.233647603(8 

Iwasaki 1.05772129(4) 0.18672220(2) 0.17428813(3) 

DBW2 3.4449465(4) -0.22085461(4) -0.10748502(5) 



TABLE XX: The coefficients ef 7) - ej' 9) (Eq. ([32])) for different actions. 
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Action eip' 1 ^ e^' 2 ^ e^ 1 ' 3 ^ 

Plaquette 0.37366536(7) -1.66446414(3) -0.5750281973(1) 

Symanzik 0.51501972(4) -1.47511786(3) -0.4769739579(4) 

TILW (8.45) 0.62806240(5) -1.34231565(2) -0.411841977(2) 

TILW (8.00) 0.64974666(4) -1.31859128(1) -0.4006364262(1) 

Iwasaki 0.82253993(3) -1.15028034(3) -0.3267471901(5) 

DBW2 1.5201736(4) -0.67309671(3) -0.1483549734(1) 

TABLE XXI: The coefficients e^' l) (Eq. ([MI)) for different actions. 

Action eip' 1 ^ e^' 2 ^ e^' 3 ^ 

Plaquette -4.05372833(7) 1.866287582(5) -0.8573692476(6) 
Symanzik -3.0228493(1) 1.59558642(1) -0.7667066321(8) 
TILW (8.45) -2.28808611(8) 1.39406610(3) -0.7029813946(7) 
TILW (8.00) -2.15494125(2) 1.357142525(7) -0.691443627(2) 
Iwasaki -1.17592792(3) 1.087913642(4) -0.6063355831(7) 
DBW2 2.0163147(3) 0.15488056(9) -0.359624204(3) 

TABLE XXII: The coefficients e^' l) (Eq. ((HI)) for different actions. 

A - • (2,1) (2,2) (2,3) 

Action e T e T e T 

Plaquette 2.3328621(2) -1.52209604(8) 0.23683195(1) 

Symanzik 2.4912319(2) -1.53694399(3) 0.26295051(2) 

TILW (8.45) 2.4578347(1) -1.49009843(6) 0.26767225(1) 

TILW (8.00) 2.44550431(4) -1.48035333(8) 0.26834751(3) 

Iwasaki 2.3441345(2) -1.4848088(1) 0.30172406(2) 

DBW2 1.3013094(2) -1.2798033(3) 0.3475909(1) 

TABLE XXIII: The coefficients e^ 2,1) - (Eq. flHD) for different actions 
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Action 


£ (2,4) 


.(2,5) (2,6) 

C rp C rp 


Plaquette 


-2.02795509(9) 


0.11808647(1) 0.07250824(3) 


Symanzik 


-1.55221265(9) 


0.04504264(4) 0.07020813(1) 


TILW (8.45) 


-1.2361662(1) 


-0.00137858(3) 0.06309676(1) 


TILW (8.00) 


-1.17754207(2) 


-0.0104623(1) 0.06169762(2) 


Iwasaki 


-0.6509124(1) 


-0.11083048(4) 0.06071268(8) 


DBW2 


1.4802111(2) 


-0.51691409(6) 0.03446776(7) 



TABLE XXIV: The coefficients 4?' 4) - ep &) (Eq. ([Ml)) for different actions. 
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Action 


.(2,7) _(2,8) 

C- rp &rp 


£ (2,9) 


Plaquette 


0.3932905(2) 1.10184617(6) 


-0.02744360(1) 


Symanzik 


0.1467998(1) 1.03762644(2) 


-0.03500219(3) 


TILW (8.45) 


-0.0094312(1) 0.97633518(7) 


-0.038311849(6) 


TILW (8.00) 


-0.0364177(1) 0.96442476(2) 


-0.03892456(3) 


Iwasaki 


-0.2049940(1) 0.88259379(5) 


-0.04896801(2) 


DBW2 


-0.3118850(6) 0.51292384(2) 


-0.07146761(1) 



TABLE XXV: The coefficients ep 7) - ep 9) (Eq. ([HI)) for different actions. 



Action 


J2.1) 

Crp/ 


-(2,4) 

Orp/ 


_(2,7) 

C-rpf 


Piaquette 


0.00047095(5) 


-0.30537822(7) 


1.4070324(1) 


Symanzik 


-0.26900984(7) 


-0.67000961(5) 


1.0574111(1) 


TILW (8.45) 


-0.31308657(6) 


-0.90858179(3) 


0.7651952(1) 


TILW (8.00) 


-0.31678921(7) 


-0.95117312(1) 


0.7107479(2) 


Iwasaki 


-0.45213470(9) 


-1.24108756(9) 


0.3465297(2) 


DBW2 


-0.8461093(2) 


-1.9354110(1) 


-0.6289363(3) 



TABLE XXVI: The coefficients efi 1} , e%f\ efi 7) (Eq. ([35])) for different actions. 



Action £^/' 8 ' ) e^/' 9 ^ 

Plaquette 0.28175492(1) 0.054718244(7) 

Symanzik 0.24663315(4) 0.06136902(2) 

TILW (8.45) 0.2319752(1) 0.06397589(2) 

TILW (8.00) 0.22947913(7) 0.064400395(7) 

Iwasaki 0.19560470(5) 0.07153771(1) 

DBW2 0.13147908(9) 0.08509400(2) 



TABLE XXVII: The coefficients ef, , ( Ec l- (S3)) for different actions. 
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